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Observation number correlation in WMAP data 
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ABSTRACT 

A remarkable similarity between the large-scale non-Gaussian pattern of cosmic microwave 
background (CMB) temperatures obtained by Wilkinson Microwave Anisotropy Probe 
(WMAP) mission and the distribution feature of observation numbers is noted. Motivated 
from such a similarity, in this work we check the WMAP data for the correlation between 
pixel temperature t and observation number N. Systematic effect of imbalance differential 
observation and significant t-N correlation in magnitude, distribution non-Gaussianity and 
north-south asymmetry are found. Our results indicate that, for precision cosmology study 
based on WMAP observations, the observation effect on released WMAP temperature maps 
has to be further carefully studied. 
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1 INTRODUCTION 

The WMAP observations provide precision data for cosmology 
study. By analyzing CMB maps from the first year WMAP 
(WMAPl) data, Tegmark et al. (2003) find both the CMB 
quadrupole and octopole having power along a particular 
spatial axis and more works (d e Oliveira-Costa et al. 2004; 
lEriksen et al. 2004j : ISchwarz et alloO^ : iJaffe et al. 20051) find 
that the axis of maximum asymmetry tends to lie close to the 
ecliptic axis. A s imilar anomaly was also found in COBE maps 
dCopi et al. 2006h . The unexplained orientation of large-scale 
patterns of CMB maps in respect to the ecli ptic frame is one of the 
biggest surprises in CMB studies (iStarkman and Schwarz 2005h . 
A notable asymmetry of temperature fluctuation power in two 
opposin g hemispheres is also found in the WM APl and COBE 
results ( lEriksen et al. 2004bl : [Hansen et al . 2004'). After the release 



tected in the WMAP3 data dAbramo et al. 20061: Ijaffe et al. 2006 ; 


CoDi et al. 20071; 


Land & Maaueiio 2007 ; Eriksen et al. 2007 ; 


Paik et al. 20071; 


Vielva et al. 20071; ISamal et al. 20081) and 



WMAP data for the correlation between pixel temperature t and ob- 
servation number A^. A systematic effect of imbalance differential 
observation and significant t-N correlation in magnitude, distribu- 
tion non-Gaussianity and north-south asymmetry are detected and 
shown in §3. We give a brief discussion on the observation effect in 
WMAP data in 55. 



2 LARGE-SCALE NON-GAUSSIAN MODULATION 

To address the large scale anomalies, such as asymmetry, align- 
ment and low / power issues detecte d in WMAP dat a with differ- 
ent techniques, the WMAP team (Spe rgel et al. 20061) describe the 
observed temperature fluctuations, t, as a Gaussian and isotropic 
random field, t, modulated by a function /(n) 

f(n) =t(n)[l + /(n)] 

where /(n) is an arbitrary modulation function. They expand /(n) 
in spherical harmonics 



WMAP5 data teemui & Hipo lito-Ricaldi 2008h as well 

These apparent anomalies, if found to be cosmological ori- 
gin, will pose a big challenge to the standard model of cosmol- 
ogy. Therefore, inspecting the effects of WMAP observation on re- 
leased data at large angular scales more carefully is worth doing. 
We show in §2 that there exists in WMAP data a remarkable sim- 
ilarity between the large-scale non-Gaussian pattern of map tem- 
peratures and the distribution feature of observation numbers. Mo- 
tivated from such a similarity, in this work we further check the 
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and use maximum likelihood technique with a Markov Chain 
Monte Carlo solver to get the best fit values of fim with Imax ~ 2 
for the WMAP3 V-band map. The top panel of Fig. [T] is obtained 
based on the best fit coefficients, showing in a unifying manner the 
large scale anomalies in WMAP temperature fluctuations which is 
the same feature that has been identified in a number of papers on 
non-Gaussianity. 

The sky coverage of WMAP mission is inhomogeneous - the 
number of observations being greatest at the ecliptic poles and the 
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Figure 1. Large-scale non-Gaussian modulation features with Imax = 2 
for the WMAP3 V-band data . Top panel: Best-fit large-scale modulation 
function /(n) for temperature map. Middle panel: large-scale feature of 
map. Bottom panel: large-scale feature of l/rj{N) map. 

ecliptic plane being most sparsely observed tainshaw et al. 2007h . 
To address the observation number distribution in the same way for 
the temperature distribution, we calculate the spherical harmonic 
coefficients with Imax = 2 for the map of with A'^ being 
number of observations per sky pixel from the WMAP3 V-band 
data. The middle panel of Fig. 1 shows the map of recon- 
structed based on the coefficients fim- Comparing the two graphs 
at the top and middle of Fig. 1, we can see that the large-scale non- 
Gaussian modulation features of WMAP temperature map and scan 
pattern being similar for each other, showing considerable connec- 
tion between pixel temperature t and observation number A'^. 

We further show in the bottom panel the reconstructed result 
for the observation fluctuation map - the map of l/a{N) where 
the rms variation o-{N) = ^ {(A'' — (N))^) calculated within a 
region of ~ 1° side dimension for each sky pixel. From this map 
we notice, unexpectedly to some extent, that observation numbers 
in released WMAP data are also fluctuated. In comparing the three 
panels in Fig.[T] the modulation pattern for the observation fluctua- 



tion map (the bottom panel) is more similar to the detected anoma- 
lies (the top panel), indicating that the fluctuation of observation 
numbers could produce additional uncertainty to the recovered tem- 
perature map. 



3 TEMPERATURE-EXPOSURE COUPLING 
3.1 t-N Correlation 

The remarkable similarity between the large-scale non-Gaussian 
pattern of WMAP CMB temperatures and the distribution feature 
of observation numbers shown in the previous section prompts us 
to check the WMAP data for the correlation between pixel temper- 
ature t and observation number A*'. To inspect the t-N coupling, for 
a sky pixel i we calculate the correlation coefficient Cj_„ (i) by 

^ E,(t.O-)-?I)(^.(j)-^) 

where the summations are taken over all WMAP pixels j (in the 
original resolution) within a spherical cap centered at the ver- 
tex i with an angular radius of 10° . We use Eq. [T] to produce 
a correlation map Cj_„(i) from the WMAP5 Q-, V- and W- 
band data separately, where the vertexes i a re defined with r5 
resolution of HEALPix pixelization scheme jGorski et al. 20051) 
in the sky sphere, if more than 20 % of pixels of a cap are in - 
side the Galactic mask KQ85 cGold et al. 20081 : iNolta et al. 2008h . 
the cap is no longer used. From a correlation map, we calcu- 
late the average of absolute correlation coefficients (|Cj_„ |) over 
the full sky, {\C^_^\)south over the South Galactic hemisphere, 
(|C't-]v \)north over the north Galactic hemisphere, and the south- 
north asymmetry ratio (|Cj_„ \) south / {\C^^,^ \)„orth- 

Now we use simulations to test the significance of the magni- 
tude oft — N correlation and its north-south asymmetry in WMAP 
data. The program synfast in HEALPix software package (avail- 
able at http://healpix.jpl.nasa.gov) can create temperature maps 
computed as realizations of random Gaussian fields on a sphere 
characterized by the user provided spherical harmonic coefficients 
of an angular power spectrum. For each studied band, we pro- 
duce 50,000 simulated temperature maps with th e synfast program 
from the best fit ACDM model power spectrum dNolta et al. 2008h 
with the beam fun ction ( Hill et al. 20081) and five-year like noise 
dLimon et al. 2003h . We compute the average absolute value of t- 
N correlation coefficient and its north-south asymmetry ratio for 
each simulated CMB map in the same way as for the WMAP5 
data. Finally, we calculate the average and its standard deviation 
from 50,000 {\C,_^\) and 50,000 {\C,_„\) south /{\C,_„\)„orth 
respectively. 

The obtained results are summarized in Table[T] From Table[T] 
we can estimate the significance of the average absolute magnitude 
oft~N correlation being 4.68cr, 4.23cr and 4.78a for Q-, V- and 
W-band respectively. Therefore, the t — N coupling in WMAP5 
data is much stronger than what expected from the ACDM model. 
In other words, the inhomogeneity of WMAP exposure may pro- 
duce notable distortion in temperatures observed for the ACDM 
CMB. Especially, for each of the three bands, none of the 50,000 
estimators (|Cj_„ |) obtained from simulation exceeds the WMAP 
values in Table. 1. This is well consistent with the ~ 4a signifi- 
cance estimation for the three bands. From Table[T]we can also see 
a considerable north-south asymmetry existing in t — A' correlation 
with significance of 2.58cr, 2.18cr and 2.21a for Q-, V- and W-band 
respectively. 
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Table 1. Average t-N correlation coefficients 







(|C^*-„I> 


<l'^t-JV \) south 
<l'^t-JV \)noTth 


Q-band 


WMAP5 
Expectation 


0.116 

0.0837 ± 0.0069 


1.25 

1.010 ±0.093 


V-band 


WMAP5 
Expectation 


0.0977 
0.0719 ±0.0061 


1.20 
0.993 ± 0.095 


W-band 


WMAP5 
Expectation 


0.0942 
0.0660 ± 0.0059 


1.21 

0.995 ± 0.097 



Furthermore, we test the Gaussianity of obtained WMAP t-N 
correlation distribution. Fig. |2] shows histograms of t-N correla- 
tion coefficients from WMAP5 Q-band data (solid line) and sim- 
ulations (crosses and dotted line) for the Galactic north and south 
hemispheres separately. From Fig. [2] we see that, for both north- 
em and southern hemispheres, the simulated t-N correlation co- 
efficients are normally distributed, and the correlation coefficients 
measured from WMAP5 data visibly deviate from the normal dis- 
tribution. The histogram of t-N correlation for the Galactic south 
hemisphere from WMAP5 Q-band data is much higher than the 
error bars at both tails, indicates that there are far more caps hav- 
ing strong t-N correlation than expected. On the contrary, the his- 
togram for the north hemisphere is rather close to the simulation 
result. These results on non-Gaussianity and north-south asymme- 
try for t-N correlation distribution are consistent with that for t-N 
correlation magnitude. 

The Kolmogorov-Smimov test is also used to compare the 
measured {F) and simulated (_Fo) cumulative distribution functions 
(CDFs). We find the maximum absolute differences (Dn) between 
the measured and simulated CDFs, Z)„ = maa;(|F — Fo|), is 0.105 
for the north hemisphere and 0. 1 85 for the south hemisphere, re- 
spectively. The threshold of Z)„ for null hypothesis F = Fo at 
a — 0.01 significance level and sample size, the used sky pixels 



N„ 



4400, is known to be 1.63/^ 



0.024. The measured 



D„ are much higher than the threshold. Therefore, the null hypoth- 
esis F = Fo has been rejected at at a = 0.01 or even higher signif- 
icance level, indicating that the measured t-N correlation deviates 
significantly from simulation. 



3.2 Effect of Imbalance Differential Observation 

We find that the above results from the WMAP5 foreground- 
cleaned Q-band data with the KQ85 mask are almost the same with 
what obtained from the data without any mask; therefore, the de- 
tected t-N correlation can not be explained by the foreground ef- 
fect. We will show that the instrumentation and observation effect 
can contribute to the observed t-N correlation and then distort ob- 
served temperature maps. 
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Figure 2. Histograms of t-N correlation coefficients. Solid line: from the 
WMAP5 Q-band data. Dotted line: expectation and la en'or bar from 50000 
simulations for the ACDM CMB. Upper panel: north hemisphere. Lower 
panel: south hemisphere. 



3. 2. 1 Differential observation 

The COBE and WMAP missions measure temperature differ- 
ences between sky points using differential radiometers consist- 
ing of plus-horn and minus-horn with a fixe d separation angle 
Oteam jSmoot et al. 19901 : iBennett et al. 20033) . The beam separa- 
tion angle of WMAP radiometers is Obeam ^ 141°. Let denote ti 
the temperature anisotropy at a sky pixel i. The raw data in a certain 



band is a set of temperature differences d between pixels in the sky. 
From A'^ observations we have the following observation equations 

^1+ — ti- = di 

t2+ — t2- = d2 

or in matrix notation 
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At = d 



(2) 



Where the scan matrix A= {a{k,i)), k — I,- - ■ ,N and i = 
1, ■ • • , L with L being the total number of sky map pixels. The 
most of elements a{k,i) = except for a{k,i — k'^) — 1 and 
a{k, i = k^) — —1, where fc^ denotes the pixel observed by the 
plus-hom and fc~ the pixel observed by the minus-horn at an ob- 
servation k. 

The normal equation of Eq.[2]is 



Mt = AM 

with M = A'^A. 

The Eq.[3]can be expressed as 



(3) 



k+-- 



(j = l,2,...,i). 

Where X]fc+=i means summing over A^^^ observations while the 
pixel i is observed by the plus-horn and X]fc-=i means summing 
over N[' observations while the pixel i is observed by the minus- 
horn, and the total number of observations for the pixel i is Ni = 
+ N~ . From the above equations we can derive the following 
iterative formula 



(n + l) 



(i = l,2,...,L). (4) 

With Eq. |4] when the number n of iteration is large enough, 
we get the final solution ii = t^"' for each pixel i. The Eq. |4] 
and the approximate iterative formula used by the WMAP team 
jHinshaw et al. 2003h . both have good performance for the differ- 
ential data of a noiseless instrument. 



3.2.2 Temperature distortion by instrument and observation 
imbalances 

The differential data of WMAP contain errors caused by the in- 
strument imbalance: the output of a WMAP radiometer from the 
plus-hom a nd mines-horn is not a purely differen tial response to 
sky signals jjarosik et al. 20031 : Ijarosik et al. 2007h . Instead of the 
J. - , a real observed differential data is 



ideal case dfc = tf,+ — t 
dk = tf.+ - t^- + Sk , 
where the difference distortion 

Sk = x^„,{tf.+ + t^-) 



(5) 



(6) 



with the transmission imbalance factor X jm being valued be tween 
about 0.001 and 0.02 for different bands jjarosik et al. 20071) . 

To roughly estimate the magnitude of the transmission imbal- 
ance effect on recovered temperatures, we use the first approxima- 
tion derived from Eq.|4]with the initials t^^^ = 

fc + — 1 k —i 

Substituting Eq.[5]into Eq.|7]we get 



(8) 



where j £ Ri, Ri is the scan-ring with angular radius 141° to the 
pixel i. For A'^^^^ observations j = 1, • ■ ■ , A'^^^, the plus-hom points 




Difference of observation number N*-N- 

Figure 3. Average temperature vs. observation number difference A'^+ — 
N~ from WMAP5 data for Q-band, V-band, W-band and ILC maps. EiTor 
bars are only marked on the graph for ILC map, which are similar for other 
graphs. 



to the pixel i and the mines-horn scans along the ring Ri; and for 
N~ observations j = A^^^ + 1, • • • , A''^, the mines-horn points to 
the pixel i and the plus-horn scans along the ring Ri. From Eq. [8] 
we can estimate the temperature distortion 



■Sr — tit , 



(9) 



where tn = tj/^i' Sr = J^jli Sj/Ni and j € Ri 
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Figure 4. Average temperature vs. observation number difference A'^+ — 
N~ from WMAP5 W-band map after differential imbalance compensation 
with a compensate factor 5 = 50 /xK. 



From Eq.|9]we can see that a hot source contained on the scan- 
ring Ri will let the ring average temperature t j? and the recov- 
ered temperature t^^' ^ 0, in other words, a hot foreground source 
might systematically make the recovered temperatures on its scan- 
ring lower. We have indeed found such systematic distortion in re- 
leased WMAP5 maps: scan-rings of hot sources are significantly 
cooled and strongest anti-correlations between pixel temperature 
and temperature o f scan-r ing with different separation angle 6 ap- 
pear at e ~ 141° jLiu & Li 200 9). 

Eq. [9] indicates that there might exist another kind of system- 
atic distortion in recovered temperatures caused by transmission 
imbalance of radiometers through observation imbalance 7^ 
A'^" . Eq. [6] shows that the difference distortion 5k of an observa- 
tion k dependents on the measured temperatures tf.+ + tf,-- Al- 
though the transmission imbalance factor Xim is rather small, but 
the difference distortion 5k can be considerable if any one of the 
measured temperatures is high enough. From Eq. |9] we see that 
for balance observations - the observation numbers of two horns 

= N~ , the difference distortion can not affect the recovered 
temperature. If the scan-ring Ri of pixel i contains hot sources and 
the observations are imbalance A''^:^ 7^ A^~, the recovered tem- 
perature might be distorted. To check this effect, we arrange the 
temperatures ti of all pixels i out of the foreground mask KQ75 
from a released WMAP5 map in order of the observation number 
difference ANi = iV + - iV " from the WMAP5 TOD data. Fig. [3] 
shows the dependence of temperature averaged over a difference 
interval of 10 vs. difference AA'^ between observation numbers by 
plus-horn and mines-horn. For Q-, V-, W-band and ILC maps, tem- 
perature distortions systematically changing along with AA'^ as ex- 
pected from above discussion on Eq. |9]are evidently exhibited in 
Fig. [3] 

To roughly estimate the average magnitude of difference dis- 
tortion, we try to correct the error by introducing a constant com- 
pensate factor 5 for the WMAP5 differential data dk to generate 



dk — dk 



(10) 



and reconstruct a new temperature map from the correlated TOD 
data d*. Fig. |4] shows the A'^^ — A^^ dependence from the new 
WMAP5 W-band map with 5 — 50 /iK, where the observation 
imbalance effect is much mitigated. Fig. [3] and Fig. |4] show that 
temperature distortions caused by the transmission imbalance of 
radiometers are remarkable and scan-rings for pixels with strong 




Correlation coefficient 
Figure 5. Histogram of AN-N coiTelation coefficients. 



imbalance observations (AA'' > 50) may generally contain hottest 
sources. 



3.2.3 N-AN correlation 

The dependence of average temperature (i) vs. observation number 
difference A A' shown in Fig.|3]indicates that the input transmission 
imbalance of WMAP differential radiometers and unequal horn 
coverage (AA'' 7^ 0) may generate the t-N correlation observed 
in released WMAP maps, if a considerable N-AN correlation ex- 
ists. To check it, we calculate correlation coefficients Cjy_^„(i) 
from WMAP5 Wl-band data with the same procedure we use to 
calculate Cj_„(i) in §3.1: the correlation calculation summations 
are taken over all WMAP pixels j within a spherical cap centered 
at the vertex i with an angular radius of 10°, the vertexes i are de- 
fined with r 5 resolution of HEALPix pixelization scheme in the sky 
sphere, if more than 20% of pixels of a cap are inside the Galactic 
mask KQ85 the cap is no longer used. Fig. [5] shows the histogram 
of correlation coefficients Cjy_^„ (i). From Fig. [5] we can see that 
strong correlation between the observation number Ni and differ- 
ence ANi = Nf' — N~ of a sky pixel i exists in observations used 
for making WMAP temperature maps. For most pixels i, the ob- 
servation number A'^; is positively correlated with the observation 
imbalance AA'^i, then the temperature distortion by instrument and 
observation imbalances demonstrated in this section should con- 
tribute to the detected t — N correlation presented in §3.1 through 
the N-AN correlation. 



4 DISCUSSION 

Unexplained large scale anomalies, e.g. the orientation of large- 
scale patterns in respect to the ecliptic frame, the north-south 
asymmetry of temperature fluctuation power etc, have detected in 
WMAP data. Their origin has long puzzled scientists. Several rad- 
ical explanations with an anisotropic cosmology have already been 
proposed. However, it is, in particular, hard to imagine that there 
exists a special direction related to the ecliptic or Galactic plane 
in the early universe. In this work, we show that the large-scale 
features of the WMAP observation number map and observation 
number fluctuation distribution are very similar to the large-scale 
non-Gaussian modulation feature in released WMAP temperature 
map (§2), and find that significant correlation between pixel tem- 
perature and observation number exists in WMAP data (§3.1), such 
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correlation can be produced by the systematic effect of WMAP in- 
strument and observation imbalances on temperature maps (§3.2). 
The detected observation number correlation should contribute, at 
least partly, to the apparent similarity between the large-scale fea- 
tures of the WMAP observation exposure and temperature maps, 
and hopefully provide a natural way to explain some large-scale 
anomalies in released WMAP temperature maps. 

Inhomogeneity of observation numbers used in WMAP map- 
making is emerged at different sky scales, which should generate 
systematic errors in WMAP temperature maps in a wide range of 
angular scale through the significant t-N correlation revealed in 
this work. To limit systematic artifacts, a large amount works have 
been performed by the WMAP team. Due to the differential nature 
of WMAP observations, it is a difficult task. As an example, though 
the effect of input transmission imbalance from radiometer nonide- 
alities has been noti ced, calib rated and mo dified by the WMAP 
team djarosik et al. 2 003: Jaros ik et al. 20071 ). we still find remark- 
ably systematic dependence of temperature vs. observation num- 
ber difference between the two horns residual in released WMAP 
maps. As shown in §3.22, hotspots in the sky can distort recovered 
temperatures of pixels on large part of the sky with a complicated 
way. It has to be pointed out that using Eq.[TO]with a constant com- 
pensation factor to correct the effect of horn imbalance is just to 
roughly estimate the average magnitude of differential imbalance 
distortion. What shown in Fig.|4]is only on the meaning of the av- 
erage. More works have to be done to find a proper approach to 
modify the effect of imbalance differential observation to recover a 
corrected temperature map. 

The real accuracy of cosmology parameters is the most impor- 
tant issue for high precision cosmology. Systematical temperature 
errors and structured noise fluctuations existed in CMB maps will 
certainly distort the angular power spectrum and the best-fit cos- 
mology parameters as well. It is obviously needed to further study 
the errors in WMAP temperature and noise fluctuation maps caused 
by the observation inhomogeneity and imbalance. The systematic 
distortions detected by us in released WMAP maps come from the 
WMAP's differential nature. The next CMB mission Planck is de- 
signed to measure the CMB anisotropy with completely different 
mode and expected to be unaffected by such kind of distortions. 



ACKNOWLEDGMENTS 

The referees are thanked for their helpful comments and sugges- 
tions and Prof. S.N. Zhang for suggestion on presentation of de- 
tected non-Gaussianity of t — TV correlation. This work is supported 
by the National Natural Science Foundation of China (Grant No. 
10533020), the National Basic Research Program of China (Grant 
No. 2009CB-824800), and the Directional Research Project of the 
Chinese Academy of Sciences (Grant No. KJCX2-YW-T03). The 
data analyzed in this work are obtained through the HEASARC on- 
line service provided by the NASA/GSFC. 



Copi, C.J., Huterer, D., Schwarz, D.J. & Starkman, G.D., 2006, 
MNRAS, 367, 79 

Copi, C.J., Huterer, D., Schwarz, D.J. & Starkman, G.D., 2007, 
Phys. Rev D, 75,023507 

Cruz M., Martinez-Gonzalez E., Vielva P. & Cayon L., 2005, MN- 
RAS, 356, 29 

Cruz M., Tucci M., Martinez-Gonzalez E. & Vielva P, 2006, MN- 
RAS, 369, 57 

Cruz M., Cayon L., Martinez-Gonzalez E. & Vielva P., 2007a, 
ApJ, 655, 11 

Cruz, M., Turok, N., Vielva,P., Martinez-Gonzalez E. & Hobson 
M., Science, 2007b, 318, 1612 

de Oliveira-Costa, A., Tegmark, M., Zaldarriaga. & Hamilton A., 
2004, Phys. Rev. D, 69, 063516 

Eriksen H.K., Hansen, F.K, Banday, A.J., Grski, K.M. & Lilje, 
PB., 2004a, ApJ, 605, 14 

Eriksen H.K., Banday A.J., Gorski K.M. & Lilje, PB., 2004b, 
ApJ, 612, 633 

Eriksen H.K., Banday, A.J., Grski, K.M. Hansen, F.K & Lilje, 

PB., 2007, ApJ, 660, L81 

Gold B. et al., 2008, arXiv:080370715 

Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, 
F K., Reinecke, M., Bartelmann, M., 2005, ApJ, 622, 759 
Hansen FK., Banday A.J. & Gorski K.M., 2004, MNRAS, 354, 
641 

Hill, R.S., et al. 2008, arXiv:0803.0570 
Hinshaw, G. et al., 2003, ApJS, 148, 63 
Hinshaw, G. et al., 2007, ApJS, 170, 288 

Jaffe T.R., Banday A.J., Eriksen H.K.,Gorski K.M. & Hansen 
FK., 2005, ApJ, 629, LI 

Jaffe T.R., Banday A.J., Eriksen H.K.,Gorski K.M. & Hansen 

FK., 2006, A&A, 460, 393 

Jarosik, N. et al. 2003, ApJS, 148, 29 

Jarosik, N. et al. 2007, ApJS, 170, 263 

Land K. & Magueijo J., 2007, MNRAS, 378, 153 

Limon, M. et al., 2003, Wilkinson Microwave Anisotropy 

Probe (WMAP) : Explanatory Supplement, Greenbelt, 

MD: NASA/GSFC; Available in electronic form at 

'http://lambda.gsfc.nasa.gov/ 

Liu, H. & Li, TP, 2009, Sci China G-Phy Mech Astron, in press; 
arXiv:0809.4160 

Nolta, M.R. et al., 2008, arXiv:0803.0593l 

Park C.G., Park C. & Gott J.R., 2007, ApJ, 660, 959 

Samal, PK., Saha, R., Jian, R & Ralston, J.P, 2008, MNRAS, 

385, 1718 

Schwarz, D.J., Starkman, G.D., Huterer, D. & Copi, C.J., 2004, 
Phys. Rev. Lett., 93, 221301 

Smoot, G. et al. 1990, ApJ, 360, 685 

Spergel D. N. et al., 2006, arXiv:astro-ph/0603449 /I 
Starkman, G.D. & Schwarz, D.J., 2005, Sci. Am., 293(2), 48 
Vielva P., Wiaux Y., Martinez-Gonzalez E. & Vandergheynst P., 
2007, MNRAS, 381, 932 



REFERENCES 

Abramo L.R., Bemui A., Ferreira I.S., Villela T. & Wuensche 
C.A., 2006, Phys. Rev D, 74, 063560 
Bennett, C.L. et al., 2003a, ApJ, 583, 1 
Bennett, C.L. et al., 2003b, ApJS, 148, 1 

Bernui, A. & Hipolito-Ricaldi, W.S., 2008, MNRAS, 389, 1453 



